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f AUTOMATIC LUNG PARAMETER ESTIMATOR 

The present invention relates to a device for determining one or more respiratory 
parameters relating to an individual. The device may include functionality for on-line 
5 continuous data collection, automatic assessment of the timing of measurements, 
automatic assessment of the next target (oxygen saturation of arterial blood (Sp02)), 
automatic assessment of the appropriate fraction of oxygen in inspired gas (FI02) settings 
to achieve the target Sp02, automatic control of the FI02, on-line parameter estimation, 
and automatic assessment of the number of measurements required. This functionality is 
10 achieved through a novel device including ventilatory equipment, blood gas analysis 
equipment and computer hardware and software. 



Furthermore, the present invention relates to a method for determining one or more 
p respiratory parameters by means of the above-mentioned device, wherein the individual is 
ll> suffering from hypoxemia or is at risk of hypoxemia. The individual may also be a healthy 
p individual. 

m 
n 

The use of the device for examination and monitoring respiratory parameters relating to 
n humans are of particular interest, but the device may also be applied to farm animals such 
J20 as pigs, or to domestic animals such as dogs. 

w 

j»j BACKGROUND 

jS 

Oxygen enters the body with inspiration and diffuses from the lungs into the blood. 

25 Subsequently the blood circulation transports oxygen to the tissues. Disorders of oxygen 
transport from the inspired air into the blood can result in a low oxygen saturation of the 
blood. These disorders in oxygen uptake include abnormal ventilation of the lung, seen in 
for example chronic obstructive pulmonary disease; abnormal oxygen diffusion in the 
lung, seen in for example pulmonary fibrosis; and abnormal perfusion (i.e. blood flow) 

30 through the lung. Estimation of parameters describing these oxygenation problems is 
important for diagnosis, monitoring and assessing appropriate therapeutic intervention. 
This is true in a wide variety of patients, from those who are automatically ventilated and 
who often require continuous supplement of oxygen, to out-patients who only suffer from 
dyspnoe during exercise. 

35 



>yO 00/45702 



In clinical practice the clinician usually relies upon simple measurements or variable 
estimates to assess the patients oxygenation problems. These include qualitative 
estimates obtained from stethoscopy or chest X-ray. They also include more quantitative 
estimates such as arterial oxygen saturation, the alveolar-arterial oxygen pressure 
5 gradient, or estimates of the "effective shunt", a parameter which describes all 

oxygenation problems in terms of a fraction of blood which does not flow through the 
lungs (Siggaard-Andersen and Siggaard-Andersen, 1985). 

Whilst the "effective shunt" is a parameter which has been used widely in the clinical 
10 literature it cannot adequately describe the- 'clinical' picture seen in patients when the 
inspired oxygen fraction is varied. This observation is illustrated in Figure 1 where the 
"effective shunt" has been estimated for a single patient at four different inspired oxygen 
fractions, and varies from 15-25%. 

€5 In contrast to the poor clinical description of oxygenation problems, detailed experimental 

.4 techniques such as the Multiple Inert Gas Elimination Technique (MIGET) (Wagner et a/., 

p 1974) have been developed which describe the parameters of models with as many as 

£k fifty lung compartments. The parameters of these models give an accurate physiological 

M 1 picture of the patient. Whilst the MIGET has found widespread application as an 

n 

gp experimental tool its use as a routine clinical tool has been somewhat limited (Wagner et 

P a/., 1987). This is largely due to the cost and complexity of the technique. 
W 

As stated previously, "effective shunt" is insufficient to describe oxygenation problems. 
Further parameters describing the patient's oxygenation problem can be obtained from 

25 data where inspired oxygen is varied, i.e. data similar to that presented in Figure 1. This 
was first recognised by Riley et al. (1951a, 1951b) and later by King et al. (1974). These 
authors used mathematical models to divide the oxygenation problem into that due to an 
alveolar-lung capillary drop in the partial pressure of oxygen, and that due to a shunt 
problem. To estimate two parameters describing the oxygenation problem requires taking 

30 measurements of blood samples and of ventilatory variables at each inspired oxygen 
fraction. Estimating lung parameters using the data from four inspired oxygen fractions 
required four blood samples, a procedure which is still rather time consuming and in some 
environments impractical. 



^VO 00/45702 PCT/DK00/00040 

3 

More recently, development of non-invasive methods for measuring the oxygen saturation 
of the blood have lead to renewed interest in estimation of parameters describing oxygen 
transport obtained by varying FI02. Andreassen et al. (1996, 1999), Sapsford et al. 
(1995), de Gray et al. (1997) and Roe et a!. (1997), have presented the use of two 

_ 5 parameter mathematical models of oxygen transport, the oxygenation problem being 

described as shunt combined with either a diffusion abnormality (Andreassen et al. (1996, 
1999)) or due to a ventilation/perfusion ( V I Q ) mismatch (Sapsford et al. (1995), de Gray 
et al (1997), Roe et al., (1997)). These model representations have been shown to 
provide identical fits to routine blood gas and ventilatory data obtained by varying FI02 

10 (Reeset al. 1997). 

The clinical relevance of the two parameter models is illustrated in Fig. 2, where increases 

in the pulmonary shunt parameter results in a vertical depression of the FI02/ Sa02 

C . - 

curve, (V-shift) aind abnormalities in the second parameter (ventilation/perfusion (V/Q) 

|| mismatch or oxygen diffusion resistance (Rdiff)) results in a lateral displacement of the 
P FI02/ Sa02 curve. Clearly, the lateral displacement of the FI02/ Sa02 curve (H-shift) is 
pi clinically a more significant problem as it describes a situation where large changes in 

oxygen saturation can occur for only small changes in FI02. In this situation the patient is 

at increased risk of an oxygenation problem. 



W 

p The two parameter model of Sapsford et al. (1995), has been shown to fit data from 

W- normal subjects; patients before and after thoracotomy (Sapsford et al. 1995, de Gray et 

al., 1997); and patients during (Sapsford et al. 1995, Roe et al., 1997), and after (Roe et 

al., 1997) abdominal surgery. Similarly, the two-parameter model described by 
25 Andreassen et. al. has been shown to fit data from normal subject and postoperative 

cardiac patients (Andreassen, 1999) and a wide range of as yet un-published results. 

Examples of these results are shown in Fig. 3. 

In contrary to detailed experimental approaches (e.g. the MIGET), these two parameter 
30 models can be used routinely in clinical practice. In particular, these techniques may find 
application in the monitoring and choice of therapeutic treatment for patients with left- 
sided heart failure, or to assess patients risk of post-operative hypoxaemia. 



Until now, estimation of oxygenation parameters has involved manual titration of the FI02/ 
35 Sa02 curve and off-line estimation of the parameter values. This is time consuming with 
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experimental times of approximately 45 minutes, not including the time required for off line 
parameter estimation. This limits the use of the method as a clinical tool. 

DESCRIPTION OF THE INVENTION 

5 

It is an object of the present invention to provide a device for estimation of one or more 
respiratory parameters including oxygenation parameters and lung parameters relating to 
an individual in which the necessary quantities for.enabling an estimation of respiratory 
parameters are collected automatically by a computer of the device so as to provide an 
10 automated estimation of said parameters. 

It is a further object to provide a device wherein the necessary measurements at varying 
oxygen levels are obtained in an at least semi-automated manner whereby the 
p experimental time for said estimation may be reduced. By reducing the procedural time 
these techniques have potential for routine clinical use. 

o 

gSJ It is a still further object to provide a device which is adapted for assessing a possible new 
Q target of the level of oxygen in the blood circulation based on the previously obtained 
hh measurement(s). 

m 



It is a yet still further object to provide a device, which is adapted for assessing an 
appropriate change in the current level of oxygen in the inspired gas to obtain a given 
K * ar 9 e * °f ^ e ' eve ' °f oxygen in the blood circulation. 



25 The use of the device on humans is of particular interest, but the device may also be 
applied to farm animals such as pigs, or to domestic animals such as dogs. 

The device might be of value in all kind of patients in which hypoxemia occurs or may 
occur. These conditions may e.g. be selected from the group comprising left sided heart 
30 failure, adult respiratory distress syndrome, pneumonia, postoperative hypoxemia, 
pulmonary fibrosis, toxic pulmonary lymphoedema, pulmonary embolisms, chronic 
obstructive pulmonary disease and cardiac shunting. 



Thus, the present invention relates in a first aspect of the present invention to a device for 
35 determining one or more respiratory parameters relating to an individual, comprising 
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a gas flow device having means for conducting a flow of inspiratory gas from an 
inlet opening to the respiratory system of the individual and a flow of expiratory gas from . 
the respiratory system of the individual to an outlet opening, 

a gas-mixing unit for supplying a substantially homogeneous gas to the inlet 
5 opening of the gas flow device, 

first supply means for supplying a first gas to an inlet of the gas mixing unit and 
having first control means for controlling the flow of the first gas, 

second supply means for supplying a second gas having an oxygen fraction 
different to the gas supplied from the first supply means to an inlet of the gas mixing unit 
10 and having second control means for controlling the flow of the second gas, 

a computer for determining said one or more respiratory parameters, 

first detection means for detecting the level of oxygen (Sa02, Sp02, Pa02, Pp02) 
in the blood circulation of the individual and producing an output to the computer 

fl 

^ accordingly, and 

y& 

fe second detection means for detecting the level of oxygen (FI02, FE'02, FE02, 

O 

^ PI02, PE'02, PE02) in the gas flow passing into or out of the respiratory system of the 

{3 individual and producing an output to the computer accordingly, 

^ the computer being adapted for retrieving and storing at least two measurements being 

H the concurrent output produced by the first detection means and the second detection 

CI 

means within a data structure, in which the two stored outputs are mutually related, in 

P data storage means associated with the computer, the at least two measurements being 

fi 

£j conducted at respective levels of oxygen in the gas flow passing into the respiratory 
system, the computer further being adapted for determining at least one respiratory 
parameter (Rdiff, shunt, V f Q , H-shift, V-shift) being descriptive of the condition of the 
25 individual, the determination being based on the at least two measurements. 

Hence, in its broadest aspect, the invention relates to a device for determining one or 
more respiratory parameters-relating to an individual. By the term "individual" is herein 
understood an individual selected from the group comprising humans as well as farm 
30 animals, domestic animals, pet animals and animals used for experiments such as 
monkeys, rats, rabbits, etc. 



By the term. "respiratory parameters" is herein understood parameters relating to oxygen 
transport from the lungs to the blood, such as parameters related to abnormal ventilation, 
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resistance to oxygen uptake from the lungs to the lung capillary blood, and parameters 
related to shunting of venous blood to the arterial blood stream. These respiratory 
parameters may be given as absolute values or relative values as compared to a set of 
standard values and the parameters may further be normalised or generalised to obtain 
5 parameters that are comparable to similar parameters measured for other individuals, at 
least for individuals of the same species. 



Thus, the computer may further be adapted for determining at least two respiratory 
parameters (Rdiff, shunt, V I Q , H-shift, V-shift) being descriptive of the condition of the 

10 individual, and said parameter(s) (Rdiff, shunt, V I Q , H-shift, V-shift) may alternatively or 

additionally be generalised parameters being comparable to similar parameter(s) 
determined for other individuals. 

\Q In a preferred embodiment, the computer of the device is further adapted for performing a 

Jt:5 procedure at least once, the procedure comprising 

P determining, based on at least two measurements, whether additional 

gj measurements are required, 

M» asserting a possible desired target defining a desired output of the first detection 

U means > 

Q0 producing a possible control data item based on the target, and 

US 

retrieving and storing, in the data structure, additional measurement results being 



the concurrent output produced by the first detection means and the second detection 



means. The control data item produced thereby may be outputted to a human operator by 
means of an output device so that the operator can adjust the level of oxygen in the 
25 inspired gas flow. Alternatively, the control data item may be used by another part of or a 
computer program within the computer or by an external control device for automatically 
control of the means for controlling the flow to the gas-mixing unit of at least one gas. 

According to a preferred embodiment of the present invention, the second detection 
30 means are arranged for detecting the level (FI02, PI02) of oxygen in the gas flow passing 
into the respiratory system, and the device further comprises 

third detection means for detecting the level (FE'02, FE02, PE'02, PE02) of 
oxygen in the gas flow paissing out of the respiratory system and producing an output to 
the computer accordingly, and 
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fourth detection means for detecting variables (Vt, f, V ) of the gas flow passing 
the respiratory system and producing an output to the computer accordingly, said output 
being sufficient for the computer to establish the volume flow of gas passing the 
respiratory system, 

5 the computer being adapted for retrieving and storing output from the third detection 
means and the fourth detection means within the data structure relating these stored 
output mutually as well as with the output from the first detection means and the second 
detection means retrieved simultaneously. This/these measurement(s) enable(s) the 
computer to estimate or establish the oxygen consumption of the individual, either 
10 implicitly as part of the estimation of respiratory parameters, or the computer may further 
be adapted for explicitly establishing, based on said measurement(s), the oxygen 
consumption (V02) of the individual. 

a 

j3 It is advantageous for the device according to the present invention that the computer is 
Qfe adapted to determine a parameter relating to an equilibrium state of the overall oxygen 
gjj uptake or consumption of the individual based on the output of at least one of the 

P detection means, to compare said parameter with a predefined threshold value and to 

P 

\& produce a control data item accordingly if said parameter exceeds said threshold value. 
? By determining whether an equilibrium state of the individual is obtained the timing of the 
steps of the procedure can be controlled efficiently and the overall time for performing the 
W procedure may be further reduced. 

p 

^ It is also advantageous if the computer is adapted to asses the appropriate change in 
oxygen level in the inspired gas (FI02) from the current oxygen level (FI02) so as to 

25 achieve a given desired target oxygen level in the blood (Sa02, Sp02, Pa02, Pp02) and 
produce a control data item accordingly so that the oxygen level can be adjusted 
according to the data item. The actual adjustment may be performed by an operator of the 
device, in which case the data item is outputted to an output device. Alternatively and 
preferably the computer is adapted to operate the control means for controlling the flow to 

30 the gas mixing unit of at least one gas, in response to said control data item relating to the 
assessed change in oxygen level from the computer so as to change the oxygen level 
(FI02) in the inspired gas flow accordingly. The data item may instead be outputted to an 
external device, which is suitable for performing an automated control of the control 
means so as to adjust the oxygen level accordingly. 

35 
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The assessment of change in oxygen level in the inspired gas may in an embodiment of 
the invention be based on a predefined set of data representing statistical distributions of 
variables stored within data storage means associated with the computer and on said 
measurements. Details of how this may be performed are disclosed in the detailed 
5 description of the invention. Alternatively, the assessment of change in oxygen level in the 
inspired gas may be based on the rate of change of the output of at least one of the 
detection means in response to a change in oxygen, level (FI02) in the inspired gas flow. 
Typically, the oxygen level is changed stepwise or following a ramp function and the 
change over time of the oxygen level in the blood circulation or the level of oxygen in the 
10 expired gas is monitored. However, monitoring of another gas, such as C0 2> or another 
variable of the patient may additionally or alternatively be employed. 

It is preferred that one gas is atmospheric air and that another of the gasses is more or 

less pure oxygen, i.e. has an oxygen fraction higher than that of atmospheric air, 
S Preferably in the range 0.85 to 1.00. Alternatively or additionally, another gas may be 
g supplied which has an oxygen fraction below that of atmospheric air, i.e. in the range of 
g 0.00 to 0.21, preferably of 0.00 to 0.05. Thereby the oxygen level of the inspired gas may 
Q be varied not only to level above that of atmospheric air but also below that level, thus 
(** providing a wide range of possible levels for performing measurements of the individual. 
% Tne 9as having a low oxygen fraction may be supplied from a source of more or less pure 
■Q nitrogen N 2 or another suitable physiologically neutral gas, such as helium H 2 , or it may be 

re-circulated expired gas from the individual, preferably after reduction of the level of C0 2 

in the expired gas. 

25 The device should ensure by means of a security arrangement that the oxygen saturation 
in the blood circulation of the individual is in the range of 65 to 100%, preferably for 
human beings in the range of 85 to 100% to avoid the risk of damage to organs. This 
condition varies for different species of animals. 

30 The first detection means is preferably arranged for detecting a variable relating to the 
saturation level of oxygen in the arterial blood stream by means of an invasive or a non- 
invasive technique, which latter is preferred. Thus, the first detection means is in an 
advantageous embodiment a pulse oximeter. Alternatively, the level of oxygen in the 
venous blood stream may be measured by means of an invasive or a non-invasive 

35 technique, the latter again being the preferred one. 
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According to a second aspect, the present invention relates to a device for determining 
one or more respiratory parameters relating to an individual, comprising 

a gas flow device having means for conducting a flow of inspiratory gas from an 
5 inlet opening to the respiratory system of the individual and a flow of expiratory gas from 
the respiratory system of the individual to an outlet opening, 

a gas-mixing unit for supplying a substantially homogeneous gas to the inlet 
opening of the gas flow device, 

first supply means for supplying a first gas to an inlet of the gas mixing unit and 
10 having first control means for controlling the flow of the first gas, 

second supply means for supplying a second gas having an oxygen fraction 
different to the gas supplied from the first supply means to an inlet of the gas mixing unit 
and having second control means for controlling the flow of the second gas, 

63 

Jqi a computer for determining said one or more respiratory parameters, 

tfi5 first detection means for detecting the level of oxygen (Sa02, Sp02, Pa02, Pp02) 

in the blood circulation of the individual and producing an output to the computer 



accordingly, and 

H 1 second detection means for detecting the level of oxygen (FI02, FE'02, FE02, 

H PI02, PB02, PE02) in the gas flow passing into or out of the respiratory system of the 

T^O individual and producing an output to the computer accordingly, 

P the computer being adapted for retrieving and storing a first measurement being the 

concurrent output produced by the first detection means and the second detection means 
within a data structure, in which the two stored outputs are mutually related, in data 
storage means associated with the computer, the computer being further adapted for 
25 performing a procedure at least once, the procedure comprising 

determining, based on data stored within the data structure, whether additional 
measurements are required, 

asserting a possible desired target defining a desired output of the first detection 

means, 

30 producing a possible control data item based on the target, and 

retrieving and storing, in the data structure, additional measurement results being 
the concurrent output produced by the first detection means and the second detection 
means. 
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According to a third aspect, the present invention relates to a device for determining one 
or more respiratory parameters relating to an individual, comprising 

a gas flow device having means for conducting a flow of inspiratory gas from an 
inlet opening to the respiratory system of the individual and a flow of expiratory gas from 
5 the respiratory system of the individual to an outlet opening, 

a gas-mixing unit for supplying a substantially homogeneous gas to the inlet 
opening of the gas flow device, 

first supply means for supplying a first gas to an inlet of the gas mixing unit and 
having first control means for controlling the flow of the first gas, 
10 second supply means for supplying a second gas having an oxygen fraction 

different to the gas supplied from the first supply means to an inlet of the gas mixing unit 
and having second control means for controlling the flow of the second gas, 
a computer for determining said one or more respiratory parameters, 
p first detection means for detecting the level of oxygen (Sa02, Sp02, Pa02, Pp02) 

jg> in the blood circulation of the individual and producing an output to the computer 
\Q accordingly, and 

g 

M second detection means for detecting the level of oxygen (FI02, FE'02, FE02, 

h* PI02, PB02, FE02) in the gas flow passing into or out of the respiratory system of the 

5: 

individual and producing an output to the computer accordingly, 
Wb the computer being adapted for retrieving and storing at least a first measurement 

C3 being the concurrent output produced by the first detection means and the second 
Q detection means within a data structure, in which the two stored outputs are mutually 

related, in data storage means associated with the computer, the computer further being 

adapted to asses the appropriate change in oxygen level in the inspired gas (FI02) from 
25 the current oxygen level (FI02) so as to achieve a given desired target oxygen level in the 

blood (Sa02, Sp02, Pa02, Pp02) and produce a control data item accordingly. 

The second aspect as well as the third aspect of the invention is disclosed above in the 
most fundamental embodiment which according to the present invention may be 
30 combined with the additional features disclosed above with relation to the first aspect of 
the invention. 



35 



The device may be used to obtain and/or compare one or more respiratory parameters 
relating to one or more individual(s). The individual may be a healthy individual, at risk of 
suffering from hypoxemia, or suffering from hypoxemia. 
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By the term "the individual is at risk of suffering from hypoxemia" is herein understood that 
the individual has a higher/increased risk of suffering from hypoxemia compared to a 
healthy individual. The increased risk of suffering from hypoxemia may e.g. be due to a 
hereditary predisposition, a post-operative condition and/or various diseases. 

By the term "hypoxemia" is herein meant that the oxygen saturation in the biood from the 
individual is below 92%. Examples of diseases that can cause hypoxemia are left sided 
heart failure, adult respiratory distress syndrome, pneumonia, postoperative hypoxemia, 
pulmonary fibrosis, toxic pulmonary lymphoedema, pulmonary embolisms, chronic 
obstructive pulmonary disease and cardiac shunting. 

The present invention also relates to a computer system comprising at least one general 
purpose computer having one or more computer programs stored within data storage 
means associated therewith, the computer system being arranged for as well as being 
adapted for determining one or more respiratory parameters according to the devices 
and/or methods disclosed above. 

Furthermore, the present invention relates to a computer program product being adapted 
to enable a computer system comprising at least one general purpose computer having 
data storage means associated therewith and being arranged suitably to determine one or 
more respiratory parameters according to the devices and/or methods disclosed above. 



GLOSSARY 



F102 



Fraction of oxygen in inspired gas. 
Pressure of oxygen in inspired gas. 

Oxygen saturation of arterial blood, measured from a blood sample. 
Pressure of ox'ygen in arterial blood, measured from a blood sample. 
Oxygen saturation of arterial blood, measured transcutaneously. 
Pressure of oxygen in arterial blood, measured transcutaneously. 



PI02 



Sa02 



Pa02 



Sp02 
Pp02 



FEC02 



Fraction of carbon dioxide in the mixed expired gas. 
Fraction of oxygen in expired gas at the end of expiration. 



FE'Q2 



FE02 



Fraction of oxygen in the mixed expired gas. 



WO 00/45702 



12 



PEC02 Pressure of oxygen in the mixed expired gas. 

PE*02 Pressure of oxygen in expired gas at the end of expiration. 

Vt Tidal volume, i.e. volume of gas breathed per breath, 

f Respiratory frequency, i.e. number of breaths per minute. 

• 5 V02 Oxygen consumption, i.e. the amount of oxygen consumed by the tissues 

per minute. 

Vd Dead space i.e. the volume of the lung not involved in exchanging gases 

with the blood. 

shunt Respiratory parameter representing the faction of blood not involved in 

10 gas exchange. 

Rdiff Respiratory parameter representing a resistance to oxygen diffusion 

across the alveolar lung capillary membrane. 
V Ventilation. 

H 

It**? 

tfj V / Q Respiratory parameter representing the balance between ventilation and 
l» 

grJ5 perfusion in a region of the lung. 

jg V-shift Respiratory parameter representing a vertical shift in plots of FI02 

Q against Sa02 , FI02 against Sp02, FE*02 against Sa02, or FE'02 
against Sp02 . 

j»* H-shift Respiratory parameter representing a horizontal shift in plots of FI02 

||o against Sa02 , F102 against Sp02, FE*02 against Sa02, or FE'Q2 against 

0 Sp02. 



law 



BRJEF DESCRIPTION OF THE FIGURES 



25 Fig.1. Plot of the inspired oxygen fraction (FI02, x-axis) against the arterial oxygen 
saturation (Sa02, Sp02, y-axis) for 1 patient. For each data point (A-D) the "effective 
shunt" has been estimated from a single parameter shunt model (Siggard-Andersen and 
Siggaard-Andersen 1985), giving values of point A = 15%, point B = 15%, point C = 20%, 
point D = 25%. 



30 



Fig. 2. Plots of the inspired oxygen fraction (FI02, x-axis) against model predicted arterial 
oxygen saturation (Sa02, Sp02, y-axis) for 1) a normal subject with shunt = 5% and Rdiff 
= 0 kPa/(l/min) (solid line), 2) a hypothetical patient with a Rdiff or ventilation/perfusion 
disorder (dotted line), and 3) a hypothetical patient with a shunt disorder (dashed line). 
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Line A illustrates the vertical displacement of the curve (V-shift) due to a shunt disorder, 
whilst line B illustrates the horizontal displacement of the curve (H-shift) due to a 
ventilation perfusion of oxygen diffusion abnormality. 

5 Fig. 3. Plots of the inspired oxygen fraction (FI02, x-axis) against arterial oxygen 
saturation (Sa02, Sp02, y-axis). Each of the vignettes illustrates data (crosses) and 
model predicted curves fitted, to this data from: A - a normal subject (shunt = 5%, Rdiff = - 
1.5 kPa/(l/min)), B - a post-operative cardiac patient (shunt = 9.5%, Rdiff = 81.0 
kPa/(l/min)), C - a post-operative hysterectomy patient (shunt = 7%, Rdiff = 15.2 
10 kPa/(l/min)), D - a poorly compensated cardiac patient (shunt = 15%, Rdiff = 22.9 

kPa/(l/min)) f and E - a patient residing in the intensive care unit (shunt = 7%, Rdiff = 31.0 
kPa/(l/min)). 

15 Fig. 4. Experimental set-up working with nitrogen for subathmospheric oxygen levels. The 
W5 system includes: 1) A Gas Delivery Unit including gas inlets (1a, 1b), a gas mixer (1c), a 

flow or pressure gradient (1d), and equipment for the measurement arid/or setting of 
PI inspired oxygen fraction (FI02), tidal volume and respiratory frequency (1e); 2) Equipment 

for measurement of expired gases including an oxygen monitor placed so as to measure 
* end tidal oxygen fraction (2a), and/or an expiratory reservoir, used with an oxygen monitor 
£p and/or a carbon dioxide monitor to measure the fraction of gas in or leaving the expiratory 

%% reservoir (FE02, FEC02) (2b); 3) Measurement of arterial oxygen saturation (Sa02) via 
P e.g. a pulse oxymeter (Sp02); 4) Measurements of arterial or venous blood gas samples 
(optional); 5) Measurement of cardiac output (optional); 6) A computer system including 
software for automatic collection of data (6a), monitoring the steady state of the 
25 patients/subjects oxygenation (6b), a feedback controller for adjusting inspired oxygen 
fraction (6c), and estimation of gas exchange parameters. Dashed arrowed lines illustrate 
the flow of information to the computer. Dotted arrowed lines illustrated the control of the 
gas delivery unit by the computer. 

30 Fig. 5. Experimental set up using a rebreathing technique for subatmospheric oxygen 
levels. Figure 5 illustrates a modification to the set-up of Figure 4. It includes all other 
components illustrated in Figure 4, plus a carbon dioxide removal device to eliminate 
carbon dioxide from the re-inspired gases (box 7). All other points 1-6 are the same as 
Fig. 4. 

35 
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Fig. 6. Flow chart for a measurement of variables for determination of lung parameters. 
A: Begin parameter estimation if FIO2>1.00 and SpO2>0.85 

B; Continuous data recording from gas delivery unit, pulse oxymeter and expiratory gas 
measurement devices. 
5 C: Set oxygen level (FI02). 
D: Monitor 02 equilibrium. 
E: Equilibrium level. 
F: Record measurement. 
G: Sufficient number of measurements? 
10 H: Estimate new FIOZ 

I: Estimate Pulmonary Parameters. 

Fig. 7. (algorithm 1) Assessing whether another measurement is necessary and 
p determining the target Sp02 for that measurement. If current FI02 = 1.00 and Sp02 < 

1§^ 0.85% do not perform measurement. 

CO 

2q A: Is there 1 measurement of (Sp02) 1 where 0.85 < (Sp02) 1 < 0.92? 
|| B: Target Sp02: 0.85 < (Sp02) 1< 0.92 
p C: Was FI02 = 1.00 at this measurement? 

D: Patient too sick for measurement. 
2bJ E: Is there 1 measurement of (Sp02) 2 where 0.92 < (Sp02) 2 < 0.95? 
y F: Target Sp02: 0.92 < (Sp02) 2 < 0.95 
§ G: FI02 = 1.00 at this measurement? 
M> H: Target Sp02: (Sp02) 1 < Sp02< (Sp02) 2 

I: Is there 1 measurement of (Sp02) 3 where 0.95 < (Sp02) 3 < 0.98? 
25 J: Target Sp02: 0.95 < (Sp02) 3 < 0.98 

K: Was FI02 = 1.00 at this measurement? 

L: Target Sp02: (Sp02) 2 < Sp02< (Sp02) 3 

M: Set FI02 = 1.00. 

30 Fig. 8 (algorithm 2) This controller uses a mathematical model of oxygen transport with 
two parameters, shunt and either diffusion resistance or V IQ mismatch. Parameters are 
implemented as stochastic variables and as such have a probabilistic distribution. 



A: Select appropriate a priori estimates for parameters 
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The patients lung parameters are represented as stochastic variables with probability 
distributions. These parameters need to be initialised with a priori distributions. If the 
patients lung parameters have been investigated previously, or if the patient belongs to a 
well-defined population there may be well-defined a priori distributions for the patient's 
5 lung parameters. 

B: Target Sp02 = first target level 

C: Update parameter estimates with measurement data. 
10 This is a Bayesian update of the parameter estimates for the measured values. The 
output of this process being revised probability distributions for the patients* lung 
parameters. 

*j D: Is the parameter probability mass distributed within range. 

03 5 If the probability distributions for the patients' lung parameters have a very narrow 

distribution, then they are estimated with good precision, and no further FI02 settings or 
measurements are required. 



« E: Predict Sp02 (distribution) when FI02 lowered/raised by a predetermined percentage, 
j* 20 using parameter estimates. The predetermined percentage is dependent on the 
(jy conditions and the patient. The mathematical models can be used to predict the effects of 
varying FI02 giving the current estimate of the probability distributions for the patients' 
lung parameters. Predictions can be obtained in terms of the probability of a certain 
oxygen saturation of the blood. 

25 

F: Is 10 % of probability mass < target Sp02. 

If the predicted probability distribution for Sp02 is distributed evenly about the target 
Sp02 then the FI02 is selected for the next measurement. 

30 G: Set the selected F102 level. 

H: Continue the algorithm only if there are more target SpQ2 levels? 



I: Set the next target Sp02 level. 

35 
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Fig. 9 illustrates a graph of a patients parameter (A, x-axis) plotted against the probability 
that this parameter takes a certain value (P(A), y-axis). One of these graphs is used for 
each patient parameter (i.e. shunt, Rdiff and or V I Q). Before a measurement procedure 
begins an a priori distribution is obtained for each of the patient parameters from 
5 computer storage. Subsequently, these a priori estimates are updated as measured data 
presents. Typical distributions of the shunt parameter are illustrated for a normal healthy 
subject both a priori (solid line, mean shunt = 5%), and following update of the distribution 
with measured data (dashed line). 



10 Fig. 10 illustrates model predicted arterial oxygen saturation (Sa02, Sp02, y-axis) when 
varying inspired oxygen fraction (FI02, x-axis). Points A and B are measured FI02/Sp02 
values which are used to update parameter values (i.e. P(parameters | measurements)). 
The updated parameter values are then used to predict the change in Sp02 on varying 
^ FI02 (i.e. P(Sp02 I FI02)). These predictions are illustrated for two different FI02 levels 
0^5 (C and D) and are plotted as probability distributions. The appropriate FI02 level is then 
Q selected so that < x% (in this case 10%) of the probability distribution is below the target 
Sp02 level (E). 



P 

3I 



DETAILED DESCRIPTION OF THE INVENTION 



^ The following description of preferred embodiments of the invention will focus on a device 
P for automating the estimation of lung parameters. This device (Automatic Lung Parameter 
^ Estimator = ALPE) enables reduction in the time taken to obtain estimates of oxygenation 

parameters, with the total time including on-line estimation of parameters taking 10-15 
25 minutes. By reducing the procedural time these techniques have potential for routine 

clinical use. This is only possible because of the substantial novelty in the ALPE which 

may include functionality for: 

1) On-line continuous data collection 

30 2) Automatic assessment of the timing of measurements 

3) Automatic assessment of the next target Sp02 

4) Automatic assessment of the appropriate FI02 settings to achieve the target Sp02 

5) Automatic control of the FI02 

6) On-line parameter estimation 

35 7) Automatic assessment of the number of measurements required 
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This functionality is achieved through a novel apparatus including ventilatory equipment, 
blood gas analysis equipment and computer hardware and software as described below. 

Description of the Automatic Lung Parameter Estimator (ALPE): 

The Automatic Lung Parameter Estimator (ALPE) illustrated in Figure 4 may be used to 
assess oxygenation parameters in any patient, with these parameters being useful for 
diagnostic or monitoring purposes. Monitoring of patients' lung parameters is of particular 
value for those patients with ongoing treatment for example those patients artificially 
ventilated or those receiving therapies for left-sided heart failure. 

The ALPE can automatically determine the parameters of models of oxygen transport. 
These parameters are obtained from numerous measurements including the FI02/Sp02 
curve, with this curve being constructed automatically by the apparatus for Sp02 varying 
between 0.85 to 1.00. 

ALPE illustrated in Fig. 4 includes the following (numbers before paragraphs refer to the 
numbers in Figure 4): 

1) A Gas Delivery Unit - This equipment includes: Two or more gas inlets, shown 
here delivering a) oxygen or nitrogen, and b) air; c) A gas mixer capable of mixing 
two input gases to the required fraction or concentration; d) A means of delivering 
the gases to the patient/subject i.e. a flow or pressure gradient; e) Equipment for 
the measurement and/or setting of inspired oxygen fraction (FI02), tidal volume 
and respiratory frequency (or minute volume). The gas delivery unit included in the 
system can either be a stand-alone device offering only this functionality, or any 
other device, which includes this functionality such as patient ventilation devices 
(respirators) commonly used for intensive care patients. Ventilatory gases are 
delivered to and removed from the patient/subject through a face mask, mouth 
piece combined with a nose clip, laryngeal endotracheal tube etc. 

2) Measurement of expired gases - Expired gases are measured using either a) An 
oxygen monitor, placed so as to measure expiratory gases and sensitive enough 
to give measurement of the end tidal oxygen fraction (FE f 02), i.e. the fraction of 
oxygen in the expired gases at the end of an expiration, b) An expiratory reservoir, 
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placed so as to capture expiratory gases during the course of the expiration, used 
in combination with an oxygen monitor and/or a carbon dioxide monitor sensitive 

enough to measure the fraction of gas in or leaving the expiratory reservoir (FE02, 

FEC02). 

5 

3) Measurement of arterial oxygen saturation (Sa02) via e.g. a pulse oxymeter 
(Sp02). 



4) Measurements of arterial or venous blood gas samples may be taken or may be 
10 monitored continuously by invasive means and put manually into the system. 

These measurements are optional. 

5) Measurement of cardiac output may be put manually into the system. This 
measurement is optional. 

18 

P 6) A computer system including software for 

■ m 
p 

H oxymeter (or any other measure of Sp02 or Sa02). 

0 b) Monitoring the steady state of the patients/subjects oxygenation. 

c) A feedback controller, which determines whether a further measurement is 
required and automatically adjusts the inspired oxygen fraction to the most 

25 appropriate level. 

d) Estimation of gas exchange parameters from the data collected. 



a) Automatic collection of data from the gas delivery unit (FI02, Vt, f), the expired 
gas measurement devices (FEXD2, FE02, FEC02 (optional)), and the pulse 



Dashed arrowed lines on Figure 4 illustrate the flow of information to the computer. Dotted 
30 arrowed lines illustrated the control of the gas delivery unit by the computer. 

A modification to the system is also included as part of this patent (Fig. 5). For 
environments where nitrogen (N2) or another physiologically neutral gas is not available 
the oxygen content of inspired gases can be reduced lower than air (FI02air = 21%) by 
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re-breathing expired gases. In this situation, in addition to all other components illustrated 
in Figure 4 a carbon dioxide removal device is included in the system to eliminate carbon 
dioxide from the re-inspired gases (box 7 Figure 5). All other points 1-6 described above 
are the same as Figure 4. 

5 

DETAILED DESCRIPTION OF THE FLOWCHARTS 



The flowcharts are provided solely to illustrate the invention by reference to specific 
embodiments. These flowcharts and the algorithms included herein, while illustrating 
10 certain aspects of the invention, do not portray the limitations or circumscribe the scope of 
the disclosed invention. 

Fig. 6 is a flowchart illustrating the processes involved during operation of the ALPE. 

Box A: After set-up of the equipment as illustrated in Fig. 4 and 5 the parameter 
estimation procedure begins. 

Box B: As part of this process the computer continuously collects data from the other 
equipment, including FI02 and Sp02 (and/or FE'02, Vt, f, FE02, FEC02). 

Box C: An initial inspired oxygen fraction is selected (FI02) and delivered to the patient. 
This is done automatically via the computer or manually by the doctor. Initially F102 is 
usually that of air (21%) but any other value of F102 can be used as the starting point for 
the experiment. At all times the patient/subject is required to have an arterial oxygen 
saturation (Sp02) greater than or equal to 0.85. The initial FI02 may therefore be set to a 
high level so as to achieve Sp02 > 0.85. 

After setting the inspired oxygen level the patients' oxygen system will take time to 
equilibrate. This usually occurs within 2-5 minutes after the perturbation. The equilibrium 
30 of the patients oxygen system is monitored automatically by the "steady state monitor" 
software in the computer. This functionality substantially reduces the time taken to 
perform a parameter estimation and is only possible because of the apparatus. 

Box D: The assessment of equilibrium can be performed using a number of algorithms, 
35 e.g. as follows: 
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1) The arterial oxygen saturation (Sp02) remains constant within a predefined range 
oyer a predefined time period. 

5 2) The difference between the fraction of oxygen in the inspired and expired gas 
remains constant within a predefined interval over a predefined time period. 

3) The calculated oxygen consumption (V02) remains constant within a predefined 
interval for a predefined time period. 

10 

The oxygen consumption (V02) is calculated automatically by the computer from the 
continuously monitored variables using the equation V02 = f (Vt-Vd) (FI02- FE'02) 

assuming or calculating a value of Vd t or using V02 = f Vt (FI02- FE02), or any variation 
P in this equation where a combination of measurements of end tidal or mixed expired 
l|| gases are used to estimate the oxygen consumption. 



P 



Box E: When equilibrium is achieved a measurement is recorded (Box F). 



¥ Box F: This measurement includes the current values of all continuously monitored 
M variables as described previously. It can also include measurements of blood gases in 
fjj from and arterial or venous blood and a cardiac output measure obtained from equipment 
p e.g. a pulmonary catheter. The last measurements are optional. 

Box G: Following a measurement it is decided either automatically by the apparatus or 
25 manually by the clinician whether a sufficient number of measurements have been 
performed, or whether to change the inspired oxygen fraction to a new level and take a 
further measurement when equilibrium is achieved. 

Box H: It is also decided either automatically by the apparatus or manually by the clinician 
30 what level of FI02 should be selected for a new measurement (if necessary). An 

experiment consists of not less than 2 measurements at varying FI02 levels, with Sp02 in 
the range 0.85-1.00. It is important that the setting of FI02 levels achieve data points with 
Sp02 well distributed between 0.85-1.00. 
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Examples of algorithms, which can be used to implement Box G and Box H are included 
in the next section. 

Box I: After an adequate set of measurements has been taken parameters are estimated 
5 which describe the patients tung function. Parameter estimation is performed 
automatically using one or more of the following algorithms: 



10 



1 ) Graphical estimation of displacement(s) of the FI02/Sp02 curve or the FE02/ 
Sp02 curve. 

Values of inspired or expired oxygen fraction can be plotted against the arterial oxygen 
saturation (Sp02) and graphical methods used to measure the horizontal (H-shift) and 
vertical displacement (V-shift) of the data (or interpolated data) from a normal reference 



range as illustrated in Figure 2. 



8^ 



O 2) Estimation of the parameters of models of oxygen transport. 

All data collected for each of the measurements can be used with mathematical 



25 



30 



models of oxygen transport to estimate parameters describing oxygenation. 
Parameters can e,g ; be estimated describing the shunting of pulmonary blood 
(shunt) and either a resistance to oxygen diffusion or a mismatch between the 
ventilation and perfusion of the lung. 

Algorithms for Automating boxes G and H in Fig. 6: 

Numerous algorithms can be devised which enable assessment of: 

a) Whether a new measurement is required. 

b) What is the target Sp02 for this measurement. 

c) What inspired oxygen fraction (FI02) setting should be used to obtain the target 
Sp02 

These algorithms include those with complete computer automation of points a-c, and 
where points a-c are assessed using clinical judgement. 
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Two examples of these algorithms are presented here. The first includes points a and b. 
The second includes points a and c, using mathematical models of oxygen transport to 
asses the appropriate FI02 setting. 

5 It should be noted that these algorithms are only illustrations of the control system of 
ALPE and that any other algorithms which can be used to assess points a, b and c are 
included in the patent application. 



Algorithm 1:This algorithm covers points a and b above, and is illustrated in a flowchart 
10 (Fig. 7). It should be noted that if the current FI02 = 1.0 and the current Sp02 is < 0.85, 
then the patient is too ill to perform a lung assessment. 

Algorithm 2: This algorithm covers points a and c i.e. it assesses whether a measurement 
is required and estimates the appropriate FI02 setting for the next measurement given a 
^5 target Sp02. The algorithm is illustrated in the flowchart Fig. 8. This algorithm uses a 
$j mathematical model of oxygen transport with two parameters. Parameters are 
implemented as stochastic variables and as such have probability distributions as 



Q illustrated in Figure 9. 

* 

pO In box A (Figure 9) the appropriate a priori estimates are obtained for the parameter 

fl 

p distributions. If the patients lung parameters have been investigated previously, or if the 

\H 

p patient belongs to a well-defined population there may be well defined a priori distributions 

p for the patient's lung parameters. Alternatively, default parameter settings can be used. 

An example illustrating probability distributions on a parameter e.g. "shunt" or diffusion 

25 resistance "RdifF is illustrated in Figure 9. 

In box B the predefined target Sp02 level is retrieved from computer storage. 



In box C the parameters' probability distributions are updated with the measured data. 
30 This is a Bayesian update of the parameter estimates for the measured values, such that 
the probability of the parameter values given the measurements 
(P(parameters [ measurements)) can be calculated from Bayes theorem i.e. 



P(parameters | measurements) = P(measurements | parameters) P(parameters) 
35 P(measurements) 
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The output of this process being revised probability distributions for the patients' lung 
parameters updated to reflect the new information obtained from the measurements. 
These probability distributions are usually somewhat narrower than the a priori estimates 
5 as illustrated in Fig. 9. 

Box D decides whether a further measurement is required. If the updated probability 
distributions for the patients* lung parameters have a very narrow distribution, then they 
are estimated with good precision, and no further FI02 settings or measurements are 

10 required. If a further measurement is required then it is necessary to find the appropriate 
FI02 setting so as to reach the next target Sp02. This is done in several steps: first the 
mathematical models are used to predict Sp02 when the FI02 level is lowered or raised 
by a predetermined percentage. The predetermined percentage is dependent on the 

J- conditions and the patient. Sp02 is then predicted using the updated parameter estimates 

$5 and the equation: 

0 • 

(Q P(Sp02 |( FI02)) = Z P(Sp02 | FI02 , parameters) P(parameters) 

^ par am 

where P(parameters) is the current joint probability of all the parameter estimates. 

The output from this procedure is a set of probability distributions about Sp02 on varying 
fi FI02 values, as illustrated in Figure 10. Next (box F), an FI02 level is selected. The FI02 
Q level is chosen such that a small fraction (e.g. 10%) of the predicted probability mass is 

below the target Sp02 (see Figure 10). Selecting an FI02 where only a small fraction of 

the predicted Sp02 probability mass is below the target is a safety feature of this 
25 algorithm. Effectively, it means that it is unlikely that the patients Sp02 will fall below the 

target value on modification of FI02. After setting the new FI02 level the Sp02 target is 

modified and the above procedure repeated. 



WO 00/45702 PCT/DKOO/00040 

24 

REFERENCES 

Andreassen, S. t Egeberg, J., Schroter, M.P., Andersen, P.T., (1996) Estimation of 
pulmonary diffusion resistance and shunt in an oxygen status model. Comput Methods 
5 Programs Biomed, vol 51, pp 95-105. 

Andreassen, S., Rees, S.E., Kjaergaard, S., Thorgaard, P., Winter, S.M., Morgan, C.J., 
Alstrup, P., and Toft, E.(1999). Hypoxemia after coronary bypass surgery modeled by 
resistance to oxygen diffusion. Critical Care Medicine, vol 27, pp 2445-2453. 

10 

de Gray, L., Rush, E.M., Jones, J.G., (1997). A non-invasive method for evaluating the 
effect of thoracotomy on shunt and ventilation perfusion inequality. Anaesthesia, vol. 52, 
pp 630-635. 

P. 

i§ x King, T.K.C, Weber, B., Okinaka, A., Friedman, S.A., Smith, J. P., Briscoe, W.A. (1974). 
sQ Oxygen transfer in catastrophic respiratory failure. Chest, vol. 65, pp 40S-44S. 

m 

P Rees, S.E., Rutledge G.W., Andersen P.T., Andreassen, S. (1997). Are alveolar block and 
^* ventilation-perfusion mismatch distinguishable in routine clinical data. In: Proceedings of 
2P> the European society of computers in anaesthesia and intensive care conference, 

n 

p; Erlangen, Germany, September 18-19, 1997. 

£j Riley, R.L., Counard A. (1951a) Analysis of factors affecting partial pressure of oxygen 

and carbon dioxide in gas and blood of the lungs: Theory. J Applied Physiol., vol 4, pp 77- 
25 101. 

Riley, R.L., Counard A., Donald, K.W. (1951b). Analysis of factors affecting partial 
pressure of oxygen and carbon dioxide in gas and blood of the lungs: Method. J. Applied 
Physiol., vol 4, pp 102-120. " 

30 

Roe P.G., Galdelrab, R. T Sapsford., Jones, J.G. (1997). Intra-operative gas exchange and 
post-operative hypoxaemia. European Journal of Anaesthesiology, vol 14, pp 203-210. 



NVO 00/45702 



25 



PCT/DKOO/00040 



Sapsford, DJ., Jones J.G. (1995). The Pi02 vs. Sp02 diagram: a non-invasive measure 
of pulmonary oxygen exchange. European Journal of Anaesthesiology, vol 12, pp 369- 
374. 

5 Siggaard-Andersen M, Siggaard-Andersen O (1995). Oxygen status algorithm, version 3, 
with some applications, Acta Anaesthesiol Scand. Vol. 39, Supp. 107, pp 13-20. 

Wagner, P.D., Saltzman, H.A., West, J.B. (1974). Measurement of continuous 
distributions of ventilation-perfusion ratios: theory. J. Appl. Physiol. Vol 36(5): 588-599. 

10 

Wagner, P.D., Hedenstierna, G., Bylin, G. (1987). Ventilation-perfusion inequality in 
chronic asthma. Am. Rev. Respir. Dis., vol. 136, pp 605-612. 



